How to detect a genuine quantum pump effect in graphene? 
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Quantum pumping in graphene has been predicted in recent years. Till date there have been no 
experiments indicating a graphene based quantum pump. This is not uncommon as in case of other 
non-Dirac behavior showing materials it has not yet been unambiguously experimentally detected. 
The reason being that in experiments with such materials the rectification effect overshadows the 
pumped current. In this work we answer the question posed in the title by taking recourse to "strain". 
We show that the symmetries of the rectified and pumped currents towards strain reversal can effec- 
tively distinguish between the two. 



I. INTRODUCTION 

The field of quantum pumping burst into the limelight be- 
cause of an ingenious experiment performed in 1999. 
Applying time dependent magnetic fields to a quantum 
dot, the authors indicated the existence of a pumped 
current in absence of any dc voltage biasi. However, 
the pumped current was seen to be symmetric with re- 
spect to magnetic field reversal. This was a rude shock 
since pumped currents being dependent on scattering 
amplitudes unlike two terminal Landauer conductance 
which depends on probabilities should possess no sym- 
metry with respect to field reversal 2,3 . The consensus 
now is that quantum pumping if at all present is masked 
by the rectified currents in the aforesaid experiment 4 . 
Thus in this situation how to isolate a genuine quan- 
tum pump effect. The way forward is to check for mag- 
netic field symmetry of the pumped and rectified cur- 
rents. The rectified currents are symmetric with re- 
spect to magnetic field reversal while the pumped cur- 
rents do not possess any definite symmetry with re- 
spect to field reversal. This feature has been exploited 
in some works 23 which dealt with this topic in the con- 
text of non-graphene systems. What about graphene? 
A lot many works have appeared in the literature about 
graphene based adiabatic quantum pumping, among 
the earliest was Ref. [5] wherein the importance of 
evanescent modes was brought about, spin polarized 
quantum pumping 6 have also been predicted, adiabatic 
quantum pumping in graphene bilayers 7 , in graphene 
Normal-Superconductor structures 8 and with magnetic 
barriers 9 . Finally, noise in graphene quantum pump has 
been dealt with in Ref. [10] and a graphene quantum 
pump in the non-adiabatic regime has been explored in 
Ref. [Tl1|. However till date there has been no experi- 
mental demonstration of the same. Since a host of appli- 
cations of graphene based quantum pumps have been 
predicted a necessary and fool proof mechanism for its 
detection should be in place. The purpose of this work 
is to provide such a tool. 

Applying an external magnetic field to check for sym- 
metries in graphene is unwieldy, in fact this is difficult in 
any system where magnetic fields aren't a part of the 
setup abinitio. Further, an external magnetic field may 
change the magnitude and direction of pumped and rec- 



tified currents which has little or no semblance with the 
original currents in absence of magnetic fields. Thus 
because of the twin difficulty (i) external control of mag- 
netic fields at nanoscale is difficult and fraught with unin- 
tended consequences and (ii) an external magnetic field 
brought only as a check for pumping/rectification may 
lead to drastic changes in the nature of pumping or rec- 
tification. We in this letter propose to do some blue sky 
thinking and go beyond magnetic fields in graphene. We 
propose to use strain as a parameter to check on the na- 
ture of pumped/rectified currents. It has been shown that 
the Dirac band structure in graphene is fairly insulated 
upto strain induced elastic deformations of upto 20%, 
beyond which a band gap appears. In this work we are 
always below the 20% threshold. Thus strain as lower 
values can be an effective check as it doest change the 
Dirac band structure of graphene. Further, strain in- 
duced control is easy and flexible since these can be 
effectively managed by a gate voltages 16 . Further more 
as has been in vogue strain in graphene induces pseudo 
magnetic fields 17-19 and via strain new applications with 
graphene have been predicted 20-22 . Therefore strain re- 
versal would automatically imply magnetic field reversal 
and therefore as is wont rectified currents would be sym- 
metric with respect to strain reversal while pumped cur- 
rents wouldn't. This is the main message of this work. 
Now how do these phenomena of rectification or quan- 
tum pumping originate. In the experiment which was 
supposed to be the first to show quantum pumping 1 
what seems to happen is that the time dependent pa- 
rameters may through stray capacitances directly link 
up with the reservoirs and thus indirectly induce a bias 
which is the origin of rectified currents and finally what is 
quantum pumping? It refers to an unique way to trans- 
port charge without applying any voltage bias. The rec- 
tified current in a two terminal setup is given byi 3 -: 

Ireot = ^~R dX 1 dX 2 (C l — - C 2 — ) (1 ) 

2k J s dX\ dX 2 

Herein R is the resistance of circuit path and is assumed 
to be much less than the resistance of the mesoscopic 
scatterer, while C\ and C 2 are stray capacitances which 
link the gates to the reservoirs, Xi and X 2 are the mod- 
ulated gate voltages. Finally, G is the Landauer conduc- 
tance which is just the transmission probability (T) of the 
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mesoscopic scatterer. The pumped current into a spe- 
cific lead in a two terminal system, is in contrast given 



as 24 
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where I Q - ^ \ JA dX\dX 1 . In the above equation, S aj3 de- 
fines the scattering amplitude (reflection/transmission) 
of the mesoscopic sample, the periodic variation of the 
parameters X A and X 2 follows a closed path in a param- 
eter space and the pumped current depends on the en- 
closed area A in (X U X 2 ) parameter space. Initially, the 
mesoscopic sample is in equilibrium and for it to trans- 
port current one needs to simultaneously vary two sys- 
tem parameters X^t) = Xi + AX X sin(wt) and X 2 (t) - 
X 2 + AX 2 sin(wf + 6), herein AX t defines the amplitude 
of oscillation of the adiabatically modulated parameter 
and 6 is the phase difference between the modulated 
parameters. In the adiabatic quantum pumping regime 
we consider the system thus is close to equilibrium 12 . 
The essential difference between the rectified currents 
and the pumped currents are while the former is bound 
to be symmetric with respect to magnetic field rever- 
sal (via, Onsager's symmetry) since the conductance 13 
and it's derivatives enter the formula, the pumped cur- 
rents would have no definite symmetry with respect to 
magnetic field reversal 2,3 since they in turn depend on 
the complex scattering amplitudes which have no spe- 
cific dependence on field reversal unless the scattering 
system possesses some specific discrete symmetries 3 . 
The rectified currents in the adiabatic quantum pumping 
regime we consider differ from that in the non-linear dc 
bias regime. In the latter the Onsager symmetry rela- 
tions are not obeyed! 4 - while in the former (from Eq. 1) 
they are obeyed. Further rectification can also be talked 
of when a high frequency electromagnetic field is ap- 
plied to a phase coherent conductor 15 . This case also 
falls into the non-linear regime. 
In our analysis we remain in the adiabatic weak pump- 
ing regime for both rectified as well as pumped currents. 
The weak pumping regime is defined as one wherein 
the amplitude off modulation of the parameters is small, 
i.e., 6Xj <*: Xj. In the weak pumping regime the rectified 
currents reduce to 



dT 



dT 



I r ec,-I recr [C 1 —-C 2 —] 

with I° ecl = we 2 sm(6)SXi5X 2 R/4jT 2 n. T is the transmit- 
tance through the scattering region, which in case of 
a two terminal set-up is the Landauer conductance G. 
Similarly the pumped current into say the left lead, L is - 



with f pump - wesm(6)5X l 6X 2 /2n, w is the frequency of 
the applied time dependent parameter, 8 is the phase 
difference between parameters and e is the electronic 
charge. 

In a previous work 4 - I tried to look beyond these mag- 
netic field symmetry properties and provide examples 
wherein the nature or magnitudes of the pumped and 
rectified currents are exactly opposite enabling an effec- 
tive distinction between the two. However in some of the 
examples used an artificial condition of an exactly equal 
magnitude of stray capacitances was assumed. Here 
the distinguisher is strain applied to a graphene layer 
and no artificial condition of equal stray capacitances is 
required. 



II. MODEL 

Graphene is a monatomic layer of graphite with a hon- 
eycomb lattice structure 27 that can be split into two tri- 
angular sublattices A and B. The electronic proper- 
ties of graphene are effectively described by the Dirac 
equation 28 . The presence of isolated Fermi points, 
K+ and K-, in its spectrum, gives rise to two dis- 
tinctive valleys. In this work we deal with a normal- 
insulator-strain-insulator-normal (NISIN) graphene junc- 
tion. We consider a sheet of graphene on the x-y plane. 
Strain is induced by depositing graphene onto sub- 
strates with regions which can be controlably strained 
on demand 16 ' 25 ! 26 by applying a gate voltage. In Fig. [T] 
we sketch our proposed system. The strained region is 
located between < x < L, while the insulators are lo- 
cated on its left, -d < x < 0, and on its right, L < x < L+d. 
The normal graphene planes are to the left-end, x < -d, 
and to the right-end, x > L + d. 
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Figure 1 : The model system a graphene layer with a strain 
applied to a specific part. Vi and v 2 are the strengths of the 
barriers to the left and right of the strained region. In the thin 
barrier limit they reduce to phases x\ and xi- Towards the 
left and right of the structure shown two leads LL and RL (not 
shown) are placed effectively at equilibrium. 



For a quantitative analysis we describe our system by 
the Dirac equation in presence of strain that assumes 
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the form 25 



H = vf*? 1 
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T, (3) 



where £ is the excitation energy, m is the wavefunction. 
Here %, v F (set equal to unity hence forth) are the 
Planck's constant and the energy independent Fermi 
velocity for graphene, while the o-'s denote Pauli ma- 
trices that operate on the sublattices A or B. Eq. [3] is 
valid near the valleys K and K' in the Brillouin zone and 
*P = [</4(r), t/4(r), 0f,(r), 0£,(r)] f is a spinor containing 
the electron fields in each sublattice and valley. In a spe- 
cific valley the Hamiltonian is H = v F a- ■ (p - j-ft). The 
electron dynamics in a specific valley is then determined 
by the Dirac equation in presence of a gauge field ft. In 
Ref.J25| using a tight binding model a theory for uniaxial 
strain in graphene has been developed. 
The gauge field ft is because of the strain induced mod- 
ulation of the nearest neighbor hopping parameter t. We 
write t = t + 6 where 6 is the strain induced modula- 
tion. The complex space dependent vector potential 
ft = ft x - ifty is given by- ft = H„&" Kn . Here n is the 
lattice index in tight binding model. The strain induced 
modulation as shown in Fig. 1 happens in the region 
< x < L and the horizontal hopping is modulated by 
an amount 5 over a finite region of length L. Thus vec- 
tor potential ft - 59{x)9(L - x) is along y direction which 
coincides with direction of translational invariance. 
Let us consider an incident electron from the normal side 
of the junction (x < -d) with energy E. For a right moving 
electron with an incident angle <p the eigenvector and 
corresponding momentum reads 



<A + = [l,e'T 



p = (E + E F ). 



(4) 



A left moving electron is described by the substitution 

— > n — (p. 

Since translational invariance in the y-di recti on holds the 
corresponding component of momentum is conserved. 
In the insulators, -d < x < and L < x < L + d, the 
eigenvector and momentum of a right moving electron 
are given by 



tf+ = [l,e*>] 



iipiOlTip, cos <p io x e 



p e u = (E + E F - Vi), (5) 



with i - 1,2 . The trajectory of the quasi-particles in the 
insulating region are defined by the angles 0, o . These 
angles are related to the injection angles by 



sin 0,u / sin — (E + E F )/(E + E F - Vi) 



(6) 



Here, we adopt the thin barrier limit 29 defined as, 
0,o, and d -* 0, while V, -» oo, such that p tI d -*Xi- ' n tne 
strained graphene layer, (0 < x < L), the possible wave- 
functions for transmission of a right-moving electron with 
excitation energy E > read 





\ 


- 


v ~-. % 


- 


N 

1 


X 


' 1 


1 

* 

y — 
s 


0.7 


3C 1 =JC/2,% 2 =0 

. . . x 1 =t'2,x 2 =-lt/2 
Xl«n/2,Xj=JC/2 


CD 
O 

C 0.6 

eg 

3 

§0.5 

U 

0.4 


V ^7 

\ / 
\ / 

« i 

\r 
ii 

1 






i 



8/E r 



Figure 2: Strain reversal symmetry of the conductance. Pa- 
rameters are E F = 1, E F L = 2 



with q = (E + £»cos(0) + 6 and £ = sin _1 (sin(0) + 
6/(E + E F )).Jo solve the scattering problem, we match 
the wavefunctions at four interfaces: 0| A= -,/ = i//u\ x =-d, 
0i/l.v=o = ^sUo. *sU = iA2/Ia=l, and 02/L=lw = <fr\ x =L+d, 
where, starting with normal graphene at left, = 0+ + 
.$ii0_, 0// = P/0//+ + qnj/u- + miif/m. + n/0//_,j = 1,2, 
^s = Ps x i's + + qs^s- + ms^s* + ns^s-, and finally for 
normal graphene at the right, = s l2 ^+. Solving these 
equations leads to the amplitude of reflection s u , ampli- 
tude of electron transmittance s l2 . 
After this we are in position to calculate the pumped and 
rectified currents as follows: 

In our case Xi represents first thin barrier x\ and X 2 
that for second thin barrier xi- To invoke pumping in 
our proposed system we modulate the strengths (xi = 
Xw +x P sin(wO) and (xi = X20 +x P sin{wt + 9)). Herein w 
is the pumping frequency and 9 is the phase difference 
between the two modulated parameters. Thus in this 
adiabatic pumping regime the system is close to equilib- 
rium. 



III. RESULTS 

The first issue we tackle is the conductance. We inte- 
grate the transmission probability over all angles of inci- 
dence. 
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The conductance is symmetric under strain reversal as 
can be seen from FigEJ We take electron excitation en- 
ergy E - throughout and E F is as mentioned in figure. 
That implies rectified currents are symmetric too. It is 
also periodic as function of the strength of the insulating 
barrier's x'i s ( not plotted here). We also plot the con- 
ductance Fig.|3]as function of the length of strained re- 
gion, The conductance is damped for increasing lengths. 
However, one must point out periodic revivals in the con- 
ductance because of a Fabry-Perot type interference 
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Figure 3: Conductance as function of dimensionless length(L) 
of strained layer. Parameters are E F = 1,6/ E F = .5 



l p , i.e. divided by / . We see unlike the conductance the 
pumped currents Fig. S as is wont are non-symmetric 
with respect to strain reversal. Since they are uniquely 
determined by amplitudes and not probabilities. We also 
plot the pumped currents Fig.[5]as function of the length 
L of the strained region and not unlike the conductance, 
pumped currents are too damped ofcourse interspread 
with periodic revivals implying role of Fabry Perot inter- 
ference because of the two barriers on either side of the 
strained layer. A comparison of Figs. |3] and [5] shows an- 
other distinction between pumped and rectified currents. 
One can clearly see the conductance is exactly equal for 
X\ - Xi - x/2 and x\ - ~X2 - n/2 but the pumped cur- 
rents are completely different as function of strain length. 



which occurs due to the double barrier set up we have 
here. Next for the pumped currents, we adiabatically 
modulate the strengths of the two thin barriers on either 
side of the strained region, xuXi- In the weak pumping 
regime we have- 
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the plots below we plot the normalized pumped currents 



IV. CONCLUSIONS 



The experimental realization of this structure is not at 
all difficult, since double barrier structures in graphene 
have been experimentally realized for long. The only 
thing is by having two ac dependent gate voltages to 
modulate the strength of the double barrier structure. Af- 
ter this we modulate the strain applied by application of 
a gate voltage again. If this condition is realized then 
this very simple structure will be a very good identifier of 
a genuine quantum pump effect if present. 
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